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Abstract 

Purpose: To investigate the possibility of localization of intraprostatic lesions (IL) with contrast-enhanced ultrasound 
(CEUS) to support the brachytherapy treatment planning of temporary implants. 

Material and methods: Two brachytherapy treatment plans were generated for 8 patients treated with external beam 
radiotherapy and pulsed-dose rate brachytherapy boost for prostate cancer. The first and second brachytherapy treat- 
ment plan was without and with knowledge of the localization of the ILs, respectively. Pairwise comparison was per- 
formed on prostate, rectum, and urethra dose-volume parameters and total reference air kerma (TRAK)-values. 

Results: Coverage of the ILs by the 140% isodose was increased from mean 66.0-67.7% for the standard plan to mean 
92.5-95.7% for the adapted plan. The mean D90 of the ILs increased from 1.49-1.57 Gy/pulse to 1.76-1.81 Gy/pulse. Dose- 
volume parameters for the prostate, rectum, and urethra and the TRAK did not change. 

Conclusions: CEUS technique is a promising method for IL localization to aid in brachytherapy treatment planning. Dose 
coverage on the IL could be improved without any increase of dose in organs at risk. 
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Purpose 

Histologic examination of the prostate harboring a ma- 
lignancy often reveals multiple areas with malignant lesions 
[1]. These malignant lesions differ in size and morphologic 
appearance and can be divided in dominant lesions and 
smaller insignificant lesions [2,3]. If radiotherapy is used for 
treatment it is a common practice to treat the whole prostate. 
In general, the dose to the prostate is evenly distributed with 
the same dose applied to malignant lesions as to healthy 
prostate areas. Van Lin et al. proposed identifying dominant 
intraprostatic lesions (DILs) better in order to treat the en- 
tire prostate to a moderate dose and giving a boost on these 
DILs [4]. To apply a boost dose to DILs, imaging techniques 
are needed to identify these areas. Several techniques are 
reported as able to identify tumor lesions within the pro- 
state [5,6]. One of these techniques is Magnetic Resonance 
Imaging (MRI). T2-weighted imaging (T2-WI) provides 
a morphological localization of tumor lesions. T2-WI can 
be combined with functional imaging for additional in- 
formation on tumor presence. For the purpose of functional 
imaging, dynamic contrast-enhanced MRI (DCE-MRI), 
diffusion-weighted imaging (DWI), and magnetic spec- 
troscopic imaging (MRSI) can be used. There is increased 
evidence that the so-called multiparametric MRI (mpMRI), 
in which different sequences and combination of morpho- 
logical and functional MRI are used, will increase the dia- 



gnostic accuracy for detection, localization, and characte- 
rization of prostate cancer. However, at the European Con- 
sensus Meeting in December 2009 it was clear that there 
is a large disagreement between experts in the field on cri- 
teria to standardize the conduct and interpretation of 
mpMRI [7]. Although very promising, for the above-men- 
tioned reasons mpMRI still cannot be regarded useful for 
the routine clinical use [8] . Positron emission tomography 
(PET) lacks the high spatial resolution MRI has. Combina- 
tion of PET with computer tomography (CT) offers better 
depiction of tumor areas within the prostate. Promising ra- 
diotracers are carbon-11 ( n C-choline and n C-acetate) and 
fluoro-18 ( 18 F-fluorocholine and 18 F-fluoroacetate) [9]. 

Another modality to visualize tumor lesions is use of ul- 
trasound (US) imaging. Particularly contrast-enhanced 
ultrasound (CEUS) imaging has the ability to detect these 
lesions [10] . This technique visualizes perfusion within the 
prostate. Because of increased microvessel density in tumor 
areas, discrimination with healthy prostate tissue is obtained. 
US contrast is provided by small encapsulated intra- 
venously injected gas bubbles, which behave as additional 
reflectors in the blood stream. In a correlation study with 
prostatectomy specimens, Sedelaar et al. observed contrast- 
enhanced power Doppler ultrasonography able to find up 
to 79% of lesions larger than 5 mm [11]. Modern contrast 
imaging technique uses the non-linear behavior of the mi- 
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crobubbles to increase selective imaging of the bubbles for 
tumor detection [12]. Although in the early studies on CELTS 
for prostate cancer detection, low sensitivity and specifici- 
ty rates were reported, modern ultrasound techniques 
and quantitative ultrasound imaging are more promising 
[11,13-15]. The recent reports indicating high accuracy for 
prostate cancer detection enable the CELTS technique to 
be used for prostate cancer radiotherapy [16,17]. Modern 
transperineal prostate brachytherapy is in general performed 
with US guidance. Brachytherapy is characterized by a he- 
terogeneous dose distribution and it should be possible to 
treat the whole prostate gland to an elective dose with 
a boost on macroscopic tumor areas. The aim of this study 
was to investigate if CEUS techniques can support 
brachytherapy treatment planning. Since the value of the 
CEUS technique for routine clinical use has not been de- 
termined yet, this study should be considered as a proof of 
principle for future use. 

Material and methods 

Ten patients with prostate cancer who underwent ex- 
ternal beam radiotherapy with pulsed dose-rate (PDR) 
brachytherapy were included in this study. The study is reg- 
istered with www.trialregister.nl (number NTR1168). 
Details of the treatment were previously reported [18,19]. 
In brief, patients were treated for the prostate and base 
of the seminal vesicles with 3-dimensional conformal ex- 
ternal beam radiotherapy to deliver a dose of 46 Gy in dai- 
ly 2 Gy fractions. Subsequently, within 1 week a transpe- 
rineal implantation was performed with flexible catheters 
for PDR brachytherapy. The Oncentra Prostate planning sys- 
tem (Nucletron B.V., Veenendaal, The Netherlands) was 
used for intra-operative treatment planning. Catheter po- 
sitions were determined for provisional treatment planning. 
Definitive treatment planning was done with CT-scan with- 
out the rectal US probe to closely simulate the actual treat- 
ment geometry. The prescribed brachytherapy dose on the 
planning target volume (PTV) was 28.8 Gy in 1.2 Gy puls- 
es and a time of 2.0 hours between pulses. The PTV was de- 
fined as the prostate gland without margin. 

CEUS image processing and identification 
of intraprostatic lesions 

Prior to start with external beam radiotherapy patients 
underwent an examination with CEUS as an outpatient 
procedure. The sulphur hexafluoride microbubbles (Sono- 
vue™ from Bracco, Milan, Italy) were injected intraveno- 
usly via the antecubital fossa. US imaging was performed 
using an iU22 US scanner (Philips Healthcare, Bothel, USA) 
and a C8-4v endocavity probe. The power modulation tech- 
nique was used for CEUS imaging. All imaging sequences 
were stored in the DICOM (Digital Imaging and Commu- 
nication in Medicine) format, and transferred to a person- 
al computer for further analysis. Off-line, the intraprosta- 
tic lesions in the CEUS sequences were delineated. This 
contouring was done in the transverse plane on areas with 
a fast or increased enhancement as compared to the main 
enhancement of the peripheral zone. In total, 2 boluses of 
2.4 ml of contrast were injected, and the inflow of contrast 



was monitored. During up to 4 minutes after the injection 
additional planes were examined using the destruction- 
replenishment technique. In this manner, the whole prostate 
was examined for suspected lesions. Imaging was performed 
and stored in the transverse orientation with a continuous 
grayscale sweep from base to apex, together with one lon- 
gitudinal maximum cross-section. The midline of each frame 
in the transverse sweep is also located in the longitudinal 
cross-section. Using correlation techniques we were able to 
correlate each plane in the transverse sweep to one location 
in the longitudinal cross-section, and 3-dimensional (3D) re- 
construction of the prostate was made. Each CEUS record- 
ing in a transverse plane was searched for in the recorded 
sweep, and therefore the location of each CEUS recording 
could be reconstructed into the 3D data set. The delineated 
lesions were drawn in this 3D data set, and the result was used 
for the fusion with the planning system imaging data set. 

Image fusion 

The US images used for brachytherapy planning were 
acquired using a 2D side viewing US probe which position 
was controlled by a stepper. The images had a resolution 
of 0.4 x 0.4 mm 2 in-plane with 1 mm slice separation. 
The brachytherapy planning images were combined into 
a 3D volume and registered to the 3D reconstructed US scan 
as described above. This data set had 0.17 x 0.17 x 0.17 mm 3 
voxels. Besides the reconstructed US scan, a corresponding 
3D data set was available in which the CEUS detected in- 
traprostatic lesion (IL) was indicated as binary mask. All 3D 
data sets were loaded into an in-house software system for 
image registration and fusion [20]. The software was ex- 
tended with a rigid registration method suitable for US-US 
registration similar to the method proposed by Roche et al. 
[21]. First, both US images were manually aligned for 3D 
translation and rotation (6 degrees of freedom). During this 
step both US data sets were visualized in 3 orientations 
(transverse, sagital and coronal) and overlaid in a green-pur- 
ple color wash or as an interactive cut display. Then, the US 
volumes were processed by an unsharp mask filter to ex- 
tract gradient information. The final registration was per- 
formed by optimizing the correlation ratio metric between 
both processed US data sets. Both for visualization and reg- 
istration purposes, the US scans were resampled using tri- 
linear interpolation. The registration accuracy was mainly 
limited by differences in deformation between both US data 
sets because of difference in applied probe pressure. For the 
visual verification the registration accuracy was estimated 
to be better than 2 mm. After approval of the registration, 
the binary volume with the delineated ILs was overlaid on 
the brachytherapy US using the 3 D translation and rota- 
tion obtained from the US registration (Fig. 1). 

A new treatment plan was created on the fused images. 
Treatment planning was performed with the Oncentra Pro- 
state planning system. For each patient two plans were 
generated. The first plan without and the second plan with 
the ILs visualized. Dose constraints were formulated for 
the first treatment plan. The prostate volume covered by the 
reference dose (RD) should be 95% or more (VIOO-p > 95%), 
the minimal dose to the 2 ml rectal volume receiving the high- 
est dose should be 0.97 Gy/pulse or less (D2ml-r < 0.97 Gy/ 
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Fig. 1. Registered images of CEUS (A) and brachytherapy planning US (B) with delineated intraprostatic lesion 



pulse), the maximum urethral dose (Dmax-u) should be 
140% RD (1.68 Gy/pulse). The second treatment plan was 
created with the above-mentioned constraints and an ad- 
ditional constraint prescribing at least 95% of the IL volume 
covered by the 140% RD (V140-U > 95%) (Fig. 2). 

The brachytherapy treatment plans with and without 
US contrast were compared to each other. Cumulative dose- 
volume histograms (DVH) were calculated and compar- 
ed. DVH-parameters used for comparison were VIOO-p, 
V150-p, D90-p, D2ml-r, Dmax-u, and total reference air ker- 
maatlm(TRAK). 

Statistics 

Pairwise comparison of the means of the DVH-para- 
meters was done by the Wilcoxon signed rank sum test. All 
tests were two-sided, and P- values < 0.05 were considered 
significant. Statistical analysis was performed with the Pre- 



dictive Analytics Software Statistics, version 18.0 for Mac 
OS X (PASW 18.0, Chicago, IL, USA). 

Results 

One patient did not show ILs on CEUS and for another 
patient fusion of CEUS-images on the brachytherapy treat- 
ment planning images was not possible due to a large vari- 
ation in the prostate shape of the two studies. These two pa- 
tients were excluded from further analysis. The median 
prostate volume was 35.9 ml (range 16.8-50.2). Four patients 
had only 1 IL and the other 4 had 2 ILs at both sides of the 
prostate. The largest of both ILs within a prostate was cod- 
ed as IL1 and the smallest - as IL2. Six of these ILs showed 
extracapsular extension (Table 1). 

Table 2 shows that the mean V140 of the ILs was increas- 
ed from 66.0-67.7% to 92.5-95.7%. Also, the mean D90 on 




Fig. 2. Brachytherapy isodose distribution without (1A) and with (IB) knowledge of the position of an intraprostatic lesion 
in the right peripheral zone. Note the coverage of the IL by the 140% isodose 
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Table L Intraprostatic lesion (IL) characteristics 



Prostate volume 


median 35.9 ml 


range 16.8-50.2 


Number of ILs 


1 


4 




2 


4 




Extraprostatic extension 


Yes 


6 




No 


6 




Volume IL1 


median 0.25 ml 


range 0.04-0.53 


Volume IL2 


median 0.09 ml 


range 0.03-0.18 


Size IL1 


median 12.9 mm 


range 6-18.1 






range 5 9-11 0 



the ILs was increased from 1.49-1.57 Gy/pulse to 1.76- 
1.81 Gy/pulse. In 4 cases, the first treatment planning was 
performed with 12 catheters and in the other 2 cases, 14 cat- 
heters were used. In all cases the same number of catheters 
was used by repositioning the catheters corresponding to 
the location of the ILs. The adapted treatment plan did not 
lead to an alteration of the DVH-parameters of the prostate, 
rectum, and urethra (Table 3). Also, a statistically non-sig- 
nificant increase of the TRAK was observed (P = 0.09). 

Discussion 

In this pilot study we have investigated the possibility 
of introducing CEUS data for prostate brachytherapy treat- 
ment planning. To our knowledge, this is the first study in- 
vestigating this technique for brachytherapy. It appeared 
that with knowledge of the position of ILs, a better treat- 
ment plan could be generated with improved coverage of 
the ILs in the high dose regions. Better coverage of ILs was 
achieved by including the intraprostatic part as well as the 



Table 2. Intraprostatic lesion (IL) mean dose-volume 
parameters for a non-adapted and adapted treat- 
ment plan 





Non-adapted 




Adapted 


V140 IL1 (%) 


67.7 




92.7 


D90 IL1 (Gy/pulse) 


1.57 




1.81 


V140IL2(%) 


66.0 




95.7 


D90 IL2 (Gy/pulse) 


1.49 




1.76 


Table 3. Mean dose-volume parameters and TRAK 
for a non-adapted and adapted treatment plan 




Non-adapted 


Adapted 


P-value 


V100 prostate (ml) 


35.7 


35.8 


0.40 


V150 prostate (ml) 


13.7 


14.2 


0.61 


D90 prostate (Gy/pulse) 


1.36 


1.37 


0.61 


D2ml rectum (Gy/pulse) 


0.84 


0.88 


0.08 


Dmax-u (Gy/pulse) 


1.51 


1.49 


0.06 


TRAK (uGy at 1 m) 


459 


473 


0.09 



extraprostatic extension into the high-dose region. The bet- 
ter coverage of the ILs resulted also in a higher dose (D90) 
in the ILs. In this small study we could not find an alteration 
in DVH-parameters of the prostate. What is more impor- 
tant, the adapted plans did not lead to higher doses in the 
rectum and urethra. We have succeeded by just reposi- 
tioning the catheters to improve the adapted plan. Repo- 
sitioning of catheters was done on inspection of the dose dis- 
tribution in order to cover the ILs adequately. The distance 
for repositioning of a catheter was usually not more than 
5-mm compared to the original plan. We have not encoun- 
tered any need for adding more catheters in our study, but 
that may be necessary in certain situations. Considering 
repositioning of catheters one should realize that this was 
a virtual repositioning. In a real implant the catheters can 
still end up in another position with an implication on 
the DVH-parameters. However, we belief that with knowl- 
edge of localisation of ILs it will influence the positioning 
of catheters to cover them optimally. 

The TRAK was calculated to investigate changes in ex- 
posure to radioactivity. A statistically non-significant in- 
crease in TRAK was observed in the adapted plans com- 
pared to the original plans. Careful attention should be paid 
to TRAK values to be informed on exposure to radioactiv- 
ity, which cannot be found in the prescribed dose or dose 
coverage to organs. Adaptation of treatment plans can eas- 
ily lead to large differences in radiation exposure and should 
be considered in future studies. 

The question still remains how accurate the CEUS tech- 
nique is for identification of tumor lesions in order to in- 
troduce this concept with brachytherapy into the clinic. 
Halpern et al. compared the results of CEUS with 12 prosta- 
tectomy specimens [14] . In this small study they found a pos- 
itive predictive value for tumor identification of only 56% 
and a sensitivity of 42% . These numbers indicate a high false 
negative rate because tumor lesions are missed. In a simi- 
lar study, Sedelaar et al. found CEUS investigation by an ex- 
perienced investigator to have a detection rate of 61 % for 
prostate tumor lesions [11]. However, they also found a de- 
tection rate of 79% for large sized (> 5 mm) lesions. 
The implication for brachytherapy is that at least the large 
sized tumor areas, clinically significant lesions, are identi- 
fied. These high-volume tumor bearing areas probably need 
the highest dose [22]. 

A new development in the use of CEUS for prostate can- 
cer detection is the ultrasound signal analysis and its quan- 
tification. The aim is to make the analysis less subjective and 
to increase the accuracy for tumor detection. At the Eind- 
hoven University of Technology and Academic Medical 
Center, a diffusion parameter was investigated and com- 
pared to prostatectomy specimens [17]. A high sensitivity 
(81.2%) and specificity (84.6%) value was found for prostate 
cancer detection. Also the ROC area (Receiver Operating 
Characteristic area) under the curve was as high as 0.909. 
This finding is very promising and needs further investi- 
gation. Other means of detecting malignant prostatic lesions 
with high spatial resolution is the use of MRI. However, also 
with the use of MRI tumor lesions can be missed or tumor 
extension can be misdiagnosed. In a systematic review with 
meta-analysis Engelbrecht et al. found, using ROC-curves, 
a test accuracy for discriminating T2 from T3 tumors of 71 % 
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[23] . In this study, a classification of T-stage was done on 
a per-prostate level and probably the accuracy on a per-le- 
sion level would have been lower. At a European consen- 
sus meeting, there was a large disagreement between ex- 
perts in the field on how to use mpMRI for prostate cancer 
imaging [7]. Consensus was reached on 67% of items related 
to imaging parameters for tumor detection and localization, 
and on only 54% of items related to imaging interpretation 
and reporting. So, even for the more clinically used imag- 
ing technique of MRI, as compared to CELTS, further in- 
vestigation should be performed as to its validity. It is ex- 
pected that both imaging techniques will play a major role 
in the future for image guided and targeted radiotherapy 
[24] . Because in general prostate brachytherapy, implan- 
tations are performed under US guidance, identification 
of tumor lesions on US is more advantageous than on MRI. 
If MRI is used for tumor localization purposes, fusion of im- 
ages must be used if the US technique is used for implan- 
tation bringing in another source of geometric uncertain- 
ties. New developments in MRI-based implantations may 
circumvent this problem [25-28]. 

In our study we used rigid matching for contrast-en- 
hanced images and US images for treatment planning. 
Because of shape deformation we encountered some diffi- 
culties in fusion of image studies. This uncertainty may ques- 
tion the reliability of this procedure in a clinical setting. 
For this reason, patients were not actually treated with 
the adaptive plan. However, we have shown as a proof of 
principle that by identification of tumor lesions within the 
prostate brachytherapy treatment plans can be adapted for 
better lesion coverage without any increase of dose in the 
OAR. CEUS is a promising technique to aid in adaptation 
of brachytherapy treatment plans. For practical clinical us- 
age a suitable side-viewing probe should be developed for 
both contrast-enhanced images and transperineal implan- 
tations. After analysis of the images, contouring of the pro- 
state and ILs is performed for brachytherapy planning. This 
way implantation is performed without the need for image 
fusion. Without such a dedicated probe, def ormable image 
fusion techniques are absolutely needed for proper local- 
ization and positioning. When the treatment planning is 
done the implantation can be performed according to the 
plan and needle placement in the vicinity of the ILs. 

Conclusions 

Identification of ILs in the prostate with CEUS can help 
with adapting brachytherapy treatment planning for an im- 
proved coverage of the ILs without increasing the dose to 
the OARs. Further developments in image fusion and US 
hard- and software are needed before introduction into 
the clinic. 
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